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ABSTRACT: Several algorithms are presented. The standard algorithm generates all N anticliques 
of a graph G with v vertices in time 0(Nv 2 ). It can e.g. be adapted to calculate the independence 
polynomial of G, to generate all maximum cardinality anticliques, or just one maximum anticlique. 
The latter was programmed using the Mathematica 6.0 code. For a random (45, 92)-graph G a max- 
imum anticlique of size 21 was found in 1.344 sec, whereas the "hardwired" Mathematica command 
MaximumIndependentSet[G] clocked in at 155838 sec, which is five orders of magnitude slower. 



[T ! ' 1. Introduction 

Throughout G will be a finite simple graph with vertex set V = {1, 2, • • • , v}. What we call an 

anticlique X in this article is also known as independent or stable set of G, i.e. no two vertices 

in X are adjacent. The maximum cardinality of an anticlique in G is denoted by a(G), and the 

Q , number of all anticliques of G is its Fibonacci number f(G). The name derives from the fact 

c/5 ' that for n-element paths P n one clearly has 

O ' 

f(P n ) = /(P n _ x ) + /(P n _ 2 ) 

<N 

More involved recurrence relations for f(G n ) can be proven for other classes {G n : n £ N} of 
graphs [PT]. 

We first present our standard algorithm which computes f(G) for arbitrary graphs. The standard 
algorithm can be adapted to produce one maximum anticlique, or all maximum anticliques, or 
all inclusion-maximal anticliques. Although a clique of G is the same as an anticlique in the 
complementary graph G c , and thus all our results about anticliques carry over to cliques, the 
specific nature of our algorithms directly applies only to anticliques. 



Here comes the section break up. In section 2 the standard algorithm is described in detail. 
Rather than generating the N = f(G) anticliques one by one, it employs the principle of ex- 
clusion to generate "multivalued rows" r, each one of which may encode a vast number N(r) 
of anticliques. The theoretic complexity is 0(Nv 2 ). Since N(r) is readily computed and the 
number of rows r is usually small compared to N, computing N is considerably faster than gen- 
erating all anticliques by expanding the rows. Refining for each row r the arithmetic necessary 
to get N(r), one obtains the numbers of /c-element anticliques (1 < k < v) and hence the 
independence polynomial I(G,x). In particular, I(G, 1) = f(G). This is done in section 2.1. 

Section 3 is devoted to modifications of the standard algorithm. For starters, in 3.1 we fix k and 
adapt the standard algorithm so as to only produce all anticliques of cardinality > k (possibly 
none if k is too high.). If desired, one can stop as soon as the first anticlique of size > k pops 
up. In particular, provided a(G) is known beforehand (as it is in many situations), choosing 
k := a(G) — 1 yields all (or one) maximum cardinality anticliques. 

What if a{G) is not known? In 3.2 a branch and bound technique yields a(G) along with an 
a (G) -element anticlique. We compare it with the standard branch and bound method for 
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integer programming. In 3.2.1 we specialize to bipartite graphs because for them a maxi- 
mum anticlique can be found by a well known and possibly competitive 0{wy/v) algorithm 
(w = number of edges). Subsection 3.3 is devoted to numerics. For random graphs G of 
various edge density our algorithm (the one taylored to one maximum anticlique) is pitted 
against Mathematica's hardwired command MaximumIndependentSet[G], respectively against 
BipartiteMatchingAndCover [G] for bipartite graphs. We close in 3.4 with a few remarks on 
generating all inclusion-maximal anticliques. 

2. Generating all anticliques 

Let us introduce the standard algorithm, more precisely standard (a, B)-algorithm on this (5, 6)- 
graph: 



For each y G V let B(y) be the set of neighbours of y, so e.g. B(2) = {1,4}. A subset X of 
the vertex set V = {1,2, ■ ■ ■ ,v} (here v = 5) is an anticlique if and only if it satisfies these v 
conditions 

(1) a£X ^ B(a)r\X = (a € V) 

Beginning with the powerset C = 2 V we shall continuously shrink C by excluding all X which 
are not anticliques, until at the end exactly the family C C 2 V of all anticliques remains. As 
with any application of the principle of exclusion (which is discussed more formally in [Wl]), 
its efficiency hinges upon the compact represent ability of the usually huge and fast changing set 
families C. In the "worst case" a subset X C V must be represented by its characteristic 0, 1- 
vector of length v, (so X = {2,4} is (0, 1,0, 1,0)) but often that can be avoided. For instance, 
we write (2, 2, 2, 2, 2) for 2 V with the understanding that each label 2 stands for "either 1 or 
0". In view of (1), we impose upon (2,2,2,2,2) the "anti-implication" 1 — > 2,4,5, that is, we 
exclude the bad X's and retain those X G (2, 2, 2, 2, 2) that satisfy (1 G X 2, 4, 5 X). Using 
the symbolism a — > B akin of (1) we define the many-valued row 

r := (a, b, 2,6,6) 

as the family of all these X's. (Throughout the algorithm at most the five symbols 0, 1, 2, a, b, 
possibly with subscripts, will occur in any one row.) Call r an unsplit son of (2,2,2,2,2). So 
far, so good. But how to sieve all X G r that satisfy the anti-implication 2 — > 1,4? We first 
split r into the disjoint union of two candidate sons: 

n := {Xer\2eX} = (0,1,2,2,2) 

r 2 := {X£r\2£X} = (a, 0,2, b, b) 

Notice that if we force 2 G X then, because of 1 — > 2, 4, 5, we also force 1 G" X. Since 1 — > 2, 4, 5 
holds for any set X with 1 G" X, the third, fourth and fifth component of X's characteristic 
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vector are now free to be 1 or 0. This explains why n = (0, 1,2,2,2). On the other hand, if 
we force 2 ^ X, then for such X's the anti-implication 1 — > 2,4,5 is equivalent to 1 — > 4,5. 
This explains why r2 = (a, 0,2,6, b). By splitting r we managed that the new anti-implication 
2 — > 1,4 trivially holds for all X £ r2, and is easily imposed upon n by switching its fourth 
component from 2 to 0: 
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A few comments are in order. We could have continued by imposing the third anti-implication 
3 — > 4 on r's proper sons r[ and r2, then 4 — > 1,2,3,5 on all arising rows, and then 5 — > 1, 4. But 
we did it the LIFO kind of way (last in, first out), i.e. we always only processed the top row of 
our working stack, and labelled the other rows with the index of their pending anti-implications 
(PA). Thus at first both r[ and r<i had PA = 3. It just so happened that r[ satisfied 3—7-4 
already. In fact it satisfied all remaining anti-implications. Hence it was removed and constituted 
the first row of a growing output stack (shown top right). The sole remaining row of the working 
stack being subjected to the condition 3—7-4 resulted in a splitting akin to before. The top row 
had PA = 4 (which held already) and then PA = 5. Having imposed 5 — > 1,4 by splitting, the 
arising two rows are put on the output stack. Continuing that way one arrives at the complete 
output stack on the bottom right which represents a total of 

f{G) = 3 + 1 + 1 + 2 + 4 = 11 
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anticliques of G. For instance, the unique maximum anticlique X = {2, 3, 5} is a member of 
the bottom row. Notice that all singleton anticliques, as well as the empty anticlique, occur in 
various rows. We sometimes refer to a collection of many-valued rows as a context. 



Theorem: Let G be a graph with vertex set V = {1, 2, • • • , v}. 

Then the standard (a, -B)-algorithm generates the N anticliques X C V in time 0(Nv 2 ) 



Proof. We impose our anti-implications3 in the order 1 — > B(l), ■ ■ ■ ,v — > B(v), where 1, 2, • • • , v 
is any fixed ordering of V. Consider a "generic" top row r of the working stack which has t — >• B 
[B = B{t)) as pending anti-implication. In order to impose it upon r we assume by induction 
that the symbol p at position t in row r is either 0, 2 or 6 (but neither a nor 1). If we think of 
the vertices as being processed from "left to right", the situation is as follows: 

r = (. ........ ) ^p^^^ ) 

B t B 

By induction we may assume that to the right of p there is no symbol a or 1, but we also must 
not write these symbols to the right of p in the upcoming process of changing row r. 

Case 1: p = or the whole of B consists of O's. Then t — > B is satisfied by all X £ r, so r 
survives unaltered. 

Case 2: p = 2. 

Case 2.1: There is a 1 within B. Then put p = since the sets X £ r that satisfy t — >■ B are 
precisely the ones with t X. 

Case 2.2: B consists only of O's and of 2's (at least one of the latter by the Case 1). Then put 
p = a and put b on all the 2's within B. 

Case 2.3: There is no 1 and at least one a or 6 on a -B-position. Consider the row 

r = (•••,1,63,04,03,01,02,62, 2, 63,64 ,0,2^ 2,^1,0, • • •), 

B t b 

which is typical in that it illustrates all four possibilities: from a\ — > 61 only a\ is in B; from 
a 2 6262 the premise and a nonempty part (possibly all) of the anticonclusion is in B; from 
°3 ~~ > b-^bs only a proper part of the anticonclusion is in B; and from 04 — > 64 only the complete 
anticonclusion is in B. A moment's thought shows that the sets X £ r that satisfy t —> B are 
precisely the ones in the (disjoint) union of r' and r", where 

r' = (---l, 63, 04, 03, ai, a 2 , b 2 , 0, 63, 64, 0, 2, 6 2 , h, 0, •••), 

r" = (..-1, 63, 2, o 3 , 0, 0, 0, 1, 0, 0, 0, 0, 2, 2, 0, •••). 

Case 3: p = 6. 



lr The ensuing standard (a, _B)-algorithm differs considerably from the almost-namesake (a, B)-algorithm in 
[Wl] that handles implications a — » B. 
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Case 3.1: There is a 1 within B. Then put p = and leave the rest. That is, unless b constitutes 
the whole anticonclusion in its anti- implication a — > b; in which case additionally change a to 2. 

Case 3.2: There is no 1 within B. Then split r into r' and r" . Here r' has p = and is 
defined exactly as in case 3.1. As to r", it has p = 1 and O's on all U-positions. This entails 
repercussions outside B which are exactly analogous to Case 2.3. But how to change the anti- 
implication a — > b- • • b whose (say) first b = p has been turned to p = 1? If a has been in B, it 
has already been set to 0; otherwise put a = now. As to the 6's in a — > 1 b- ■ ■ b which have 
not been put to (being in B), switch them to 2. 

The above shows that one is never forced to delete rows, and that imposing an anti-implication 
on a row of length v costs 0(v). There are at most N final rows (the ones on the output stack) 
because they are mutually disjoint and each encodes at least one anticlique. Each final row 
arises by the imposement of exactly w anti-implications, and each imposing costs 0(v). Since 
no deletions occur, no other (superfluous) work has been done, and so the claimed 0(Nv 2 ) 
bound results. ■ 

Notice that the number of final rows is usually much smaller than N since one final row can 
comprise a great many of anticliques. Also, it is easy to see that the working stack can at most 
contain w rows at any given moment. If a row does not split (cases 1, 2.1, 2.2, 3.1), we say that 
it undergoes trivial changes. 



2.1. Calculating the independence polynomial 

The independence polynomial [LM1] of a graph G is defined as 



where Sk is the number of /c-element anticliques of G. Notice that f(G) = I(G, 1). As to the 
computation of the coefficients s^, each is the sum of all 



where r ranges over all finalized rows produced by the standard (a, Z?)-algorithm. Put c := 
\ones(r)\ + \twos(r)\, and let /3i, ■ ■ ■ ,(3t be the lengths of the anti-conclusions appearing in r. In 
order to get Sfc(r), first list all integer component solutions x = (x±, • • • ,x t ) of 

(2) x\ -I hx t = k-c (l<Xi<Pi, l<i<t) 

which is an easy matter. The weight of a solution x is defined as 



It is easy to see that Sfe(r) is the sum of all w(x) where x ranges over all solutions of (2). For 
instance, consider the row r below: 



I(G,x) 



a(G) 

sk%k i 



k=0 



s k (r) := \{X G r : \X\ = k}\ 



w(x) := w{x\)w(x2) ■ ■ ■ w(x t ), w{xi) : 
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Here c = 3, f3\ = 3, fa = l,/?3 = 2. Putting say k = 8 the solutions of (2) are x = (3, 1, 1) and 
x' = (2,1,2), hence 

s 8 (r) = w(x) + w(x') = Q(l + l)(l + 2)+Q(l + l)Q = 12 



It follows from the Theorem that computing the independence polynomial of a -u-graph G can 
be done in time 0(f(G)v 2 ). 



3. Variations of the standard (a, i?)-algorithm 

In order to tweak our standard algorithm in various ways we need a bit more notation. Consider 
a generic 5-valued row r and this X £ r: 

r := (0, ■■■,(), 2, •••,2, ay, 6i , a s , b s ,---,b s ) 

(3) 

X := (0,-0, 1, -..,1, 1,---,1, 0, 1,--,1, , 0, 1,---,1) 

Like symbols are adjacent in (3) only for easier visualization and notation. It would be possible 
but distracting to give a completely precise definition of a 5-valued row r = (pi, p2, ■ ■ ■ , p v ). 
Rather, the following concepts will do. Let zeros(r) be the set of positions of 0's, so 

zeros(r) := {i G V\ p\ = 0}, 

and define ones(r) and twos(r) similarly. Furthermore, let premset(r) be the set of positions 
k occupied by the symbols cij (1 < i < s). Thus premset(r) collects all "premises" a, t . Finally, 
for each k G premset(r) let anticonc(r, k) be the set of positions occupied by the symbols bi 
(assuming that at position k is aj). Formally, the X in (3) is X = V — (zeros(r) L)premset(r)), 
and it is clear that X is a maximum cardinality member of r. It is unique iff all anticonc(r, k) 
(k 6 premset(r)) have cardinality > 2. If the cardinality is 1, one can alternatively remove the 
unique (position of) bi instead of Oj. For any row r put 

(4) w max (r) := max{|X| : X G r} = v — \zeros(r)\ — \premset(r)\ 

3.1. Finding all (or one) anticliques of lower bounded cardinality 

For any natural number k we can generate the possibly empty family of all anticliques of car- 
dinality > k as follows. By induction, anchored at r = (2, 2, • • • , 2), we can assume that each 
row r that reaches the top of the working stack when running the standard algorithm, has 
w max (r) > k. If say r has PA = 12, then the 12th anti-implication is imposed upon r and gives 
rise to either an unsplit son r' or two proper sons r^,r 2 as in section 2. Say the latter happens 
(the case with r' is similar), and say i/; max (r^) < k but w ma , x {r' 2 ) > k. Then r[ is thrown away 
(since all its proper sons will also have w max < k) but r' 2 is put on top of the stack with its PA 
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set to 13. Proceeding this way the surviving finalized rows will exactly comprise the anticliques 
of cardinality > k. If we only need one such anticlique, we stop as soon as a finalized row r with 
w max (r) > k comes up. 

In particular, provided a(G) is known beforehand, one gets all maximum cardinality anticliques 
by setting k := a{G) — 1. For instance, there is a well known formula of Ryser to count all 
systems of distinct representatives. The (a, B)-algorithm can be used to generate them because 
they can be viewed as maximum anticliques in some obviously defined graph; furthermore a(G) 
is known beforehand. Another benefit of having all maximum anticliques of G concerns the core 
of G which is defined as the intersection of them all [LM2] . 

3.2 Computing one maximum cardinality anticlique 

The approach of 3.1 needs to be refined if a(G) is unknown. What we say below easily extends 
to the computation of all maximum (cardinality) anticliques and thus will not be spelled out 
(but see 3.3 for some numerics). 

To find one maximum anticlique, we initialize the output stack by (0, 0, • • • , 0) and put 
currentmax := 0. Throughout the whole procedure the output stack contains only one member. 
We change the standard algorithm as follows. Besides its PA pointer each row r in the working 
stack is accompanied by its continuously updated w mSLX (r) value. In particular, the first row 
r = (2, 2, • • • , 2) has w mscK (r) = v, and so w max (r) > currentmax. Generally, this is what 
happens when r reaches the top of the working stack. If w mSLX (r) < currentmax, then r is 
deleted. Otherwise r is processed and either gives rise to an unsplit son r', or to two candidate 
sons ri,r2- Now r' respectively r\,T2 are processed as usual but additionally each is facing 
deletion if its w max -value is < currentmax. Suppose either r' or one (or both) proper son r- 
is finalized, i.e. all anti-implications have been imposed on it. If its w nmx is < currentmax, 
then it is thrown away. On the other hand, if say w max (r 2 ) > currentmax, then we update 
currentmax := w iaax (r' 2 ) and substitute the previous member of the output stack by r 2 . The 
algorithm stops when the working stack is empty. At this moment the unique row r max of the 
output stack is such that every (often unique) maximum cardinality member X of r max is a 
maximum anticlique of G. Thus a(G) = \X\ is found. 

What we just described is a branch and bound algorithm, albeit of a novel kind with its 0, 1, 2, a, b 
symbolism. Let us compare it with, say, the standard branch and bound procedure for integer 
programming (IP). As to branching, that is cheap for IP. Namely, some fractional component, say 
x* A = 6.8, of the relaxed linear programming (LP) solution {x\, ■ ■ ■ , x*) yields two subproblems 
by adding the inequalities X4 < 6 and X4 > 7 respectively. As seen in the proof of the Theorem, 
branching (i.e. row splitting) is a bit more expensive in our algorithm but still benign. As to 
upper bounding, that is costly for (IP) where an LP-maximal solution needs to be computed, 
whereas evaluating (4) is almost gratuitous. As opposed to IP branch and bound, it does not 
make sense to pick the pending subproblems (= many-valued rows) merely according to their 
high 

w max( r ) value. Equally important is it that r has few anti-implications pending, i.e. PA(r) 
must be high. The present implementation always picks the top row in the working stack, but 
it is conceivable that selecting r according to another policy would be beneficial. 

Notice that the (a, -B)-algorithm can be adapted to find a maximum weight anticlique with 
respect to a weighing w : V — > N. Furthermore, the (a, l?)-algorithm could easily be parallized 
[XI]: Each processor sends its updated currentmax values to the shared memory which in turn 
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forwards it to all processors. 



3.2.1 The bipartite case 

In each graph G the complements V — X of the anticliques X are exactly the vertex covers Y, 
i.e. each edge of G is incident with at least one vertex of Y. It is easily seen that running 
the (a, i?)-algorithm (standard version or as in 3.2) it suffices to inflict all anti-implications 
y — > B(y) where y ranges over a vertex cover of G. In general it may be hard to find a small 
vertex cover but each bipartite graph has one of cardinality < ^\V\, namely the smaller of the 
two color classes. 

For any graph G let v{G) be the maximum size of a matching in G, and let r(G) be the minimum 
size of a vertex cover of G. Then v{G) < t{G) [S, p. 260], and v{G) = t{G) takes place whenever 
G is bipartite. A way to compute a(G) in a -u-element bipartite graph G is thus 

a(G) = v-t(G) = v-u(G) 

Here v{G) can be computed in time 0{wy/v) where w is the number of edges [S, Theorem 16.4]. 
Once a maximum matching is found one can construct from it a minimum vertex cover (and 
whence a maximum anticlique) in time 0(w) [S. Thm.16.6]. 

Our (a, l?)-algorithm of 3.2 benefits twofold from G being bipartite. Say the color classes have 
v\ and V2 vertices respectively, where v\ < V2 and v\ + V2 = v. First, the (a, i?)-algorithm can be 
run on a fi-element vertex cover (namely the smaller color class), and second we can start with 
currentmax := v% instead of currentmax := as before. A comparison of these two methods 
to compute a(G) is part of 3.2.2. 

3.3 Numerics 

We compared the (a, i?)-algorithm, implemented in the Mathematica 6.0 high level code, with 
the hardwired command MaximumlndependentSet. Mathematica 6.0 allows one to generate ran- 
dom w-graphs G of any desired edge density < d < 1. For each choice of v E {30, 40, 50, 60, 70, 
80,90} and d G {0.1,0.3,0.5,0.7,0.9} we generated three random (v , w)-graphs G\,G2,Gs of 
edge density d and picked the Gi for which the quotient of processing times was the least disas- 
trous for MaximumlndependentSet - unless the running time of the latter was more than 20000 
sec, in which case the other runs were omitted. Besides a(Gi), the running time^l and their 
proportions, we also report the number rsp of row splittings occuring in the (a, S)-algorithm. 
The number of times that rows only underwent trivial changes (as defined in section 2) was 
usually 20% to 60% of rsp. 

2 It goes without saying that in thousands of test runs both algorithms agreed upon q(G). All times over 1 sec 
are rounded to full seconds, and dito the proportion of running times is rounded to an integer. I thank Giinter 
Pilz, Helmut Prodinger, Stefan Wagner, and Yves Semegni for testifying many of these proportions. 
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(80,951) 


13 


75533 


56 


24587 


441 




(90, 1155) 


14 


200 753 
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What strikes the most, for v fixed the size of a(G) effects MaximumlndependentSet vastly more 
than the (a, f?)-algorithm. For instance, for v = 40 and a(G) = 7 respectively a(G) = 16, 
MaximumlndependentSet took 11 respectively 8058 seconds, whereas the times were 0.265 re- 
spectively 1.485 seconds for the (a, i?)-algorigthm. Put another way, for v fixed the performance 
gap gets the bigger the sparser G is. 



Let's conversely keep the density d fixed and see how the two algorithms compare when v is 
raised. For instance, a random graph G with v = 400 and d = 0.9 had a{G) = 5 (rsp = 7050) 
and the times were 19 sec versus 30'843 sec. For v = 3000 and d = 0.95 the (a, i?)-algorithm still 
bravely delivered a(G) = 6 (rsp = 309818) in 51161 sec; forget about MaximumlndependentSet. 

For v = 3000 graphs with density smaller than 0.95 would outdo the (a, S)-algorithm because rsp 
explodes. Interestingly however, one sporadic huge anticlique in such a graph (here "artificially" 
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inserted, of cardinality 100 or 500) is detected very fast: For v = 3000 and d = 0.95 it took 1285 
sec to find a(G) = 100, and only 51 sec to find a(G) = 500. 

As to the /c-threshold method of 3.1, for the (90, 1155)-graph G we set k := 13 and got a 14- 
element anticlique in less than 1 second! The catch is that we do not know whether this is the 
maximum size. Putting k := 14, the threshold-variant took 107 seconds to respond that all 
anticliques have size < k. Thus a{G) = 13 could be established in 107 + 1 = 108 seconds which 
is better than the 151 seconds of the currentmax-method of 3.2. Unfortunately, nobody tells us 
in advance that k = 13, 14 are the right values to pick. 

With the /c-threshold method finding all maximum anticliques is not much harder than finding 
one maximum anticlique. Specifically, the currentmax-method needed 0.156 sec to find a{G) = 3 
for the (70,2169)-graph in the table. Applying the fe-threshold method with k = 2 then delivered 
the total of 46 maximum anticliques in merely 0.234 sec. Similarly the time to get all 6 maximum 
anticliques (of size 21) for the (45, 92)-graph was 1.344 + 2.328 sec. 

Let us compare the (a, I?)-algorithm with the MATHEMATICA-command 

BipartiteMatchingAndCover which for a bipartite graph G computes both a maximum match- 
ing and a minimum vertex cover, as outlined in 3.2.1. We fixed one color class to be of cardinality 
v\ = 50, and the other had cardinality V2 = 50, 150,450, 1950 respectively. The edge densities 
were d = 0.1,0.5,0.9. Here are the times: 



V2 


d 


(a, 5)-algorithm 


BipartiteMatchingAndCover 


50 


0.1 


6438 


1 




0.5 


2 


4 




0.9 


0.125 


6 


150 


0.1 


7 


3 




0.5 


0.032 


11 




0.9 


0.031 


20 


450 


0.1 


0.156 


21 




0.5 


0.063 


36 




0.9 


0.063 


62 


1950 


0.1 


0.125 


127 



For our random G's we always got a(G) = V2 but this is irrelevant here. The key observation is 
that for fixed V2 increasing the edge density d impedes 

BipartiteMatchingAndCover since this algorithm depends linearly on the number of edges w. 
But it benefits the (a, -B)-algorithm because whenever an anti-implication i — > B(i) causes a 
row r to split into r\ (with at position i) and r2 (with 1 at position i), the row r2 collects a 
large number of 0's (namely |-B(i)| many) which in view of (4) likely makes w max (r2) smaller 
than the big (> V2) value of currentmax, and so T2 is cancelled. Indeed, e.g. starting with 
currentmax := instead of currentmax = i>2 increases the time from 0.063 to over 9 seconds 
for the bipartite graph with V2 = 450 and d = 0.5. The above also entails, and the table 
confirms it, that the sheer increase of t>2 (edge density being fixed) already benefits the (a, B)- 
algorithm. Of course, as testified by the 6438 versus 1 second, the flip side of the medal is that 
the (a, f?)-algorithm is severly beaten on sparse bipartitie graphs with roughly equally sized 
color classes. 
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3.4. Finding all inclusion- maximal anticliques 

How does one sieve out the iV max inclusion-maximal (dually: inclusion-minimal) sets among N 
subsets of {1, 2, • • • , n}? Comparing each of the N sets Xi with the pile of currently maximal 
sets, and accordingly adjusting that pile, costs 0(N 2 n). If N is much larger than n (which is 
the case in 3.4 and many other instances of the principle of exclusion) there is a faster way - 
though the theoretic bound still is 0(N 2 n). Namely, rather than manipulating N short sets Xi 
(of cardinality m n) one manipulates n long index sets containa] (of cardinality rj N) defined 
by 

(5) Containfa] := {1 < i < N : a G Xi} 

It is fairly obvious how this is to be done; see [W2, p. 113, 114] for more details. 

Coming back to anticliques, consider again the row r in 2.1. The inclusion-maximal sets con- 
tained in r obviously are these 2 3 sets: 

Xi = {1, 3, 4} U Ui U U 2 U U 3 , where 

C/i = {5} or {6,7,8} 

U 2 = {9} or {10} 

1/3 = {11} or {12,13} 

Using the method that employs the lists in (5) one can build up the pile of maximal anticliques 
by weaving in or deleting the row- wise maximal sets. One can conceive examples where a set X 
which is a maximal member of a row r is contained in X' which is maximal within another row 
r' . But this never occured in our tested random graphs. If we could simply pile up all row-wise 
maximal sets without worrying about (5), we would have an algorithm that produces the iV max 
maximal anticliques in time 0(N m3iX v 2 ) (by the Theorem). This would be better than the best 
theoretical bound [TTT]. 

The inclusion-maximal anticliques Xi of G for instance feature in the chromatic number of G, 
since this is the minimum number x(G) of Xj's necessary to cover the vertex set V. Because 
it doesn't harm to have a few superfluent non-maximal sets, it pays off to simply generate 
all row- wise maximal sets Xj. Computing x(G) in this way compares very favorably to the 
hardwired Mathematical command ChromaticNumber[G], despite the nasty bottleneck of finding 
a minimum cover by these X^s. 
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